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a Abstract. SQUID magnctomctry and polarised neutron reflectivity measure- 
ments have been employed to characterise the magnetic properties of U/Fe, U/Co 
and U/Gd multilayers. The field dependence of the magnetisation was measured 
' O ' at 10K in magnetic fields from -70kOe to 70kOe. A temperature dependent study 

Ch ' of the magnetisation evolution was undertaken for a selection of U/Gd samples. 

O , PNR was carried out in a field of 4.4kOe for U/Fe and U/Co samples (at room 

O i temperature) and for U/Gd samples (at 10K). Magnetic 'dead' layers of ~ 15A 

were observed for U/Fe and U/Co samples, consistent with a picture of interdif- 
fused interfaces. A large reduction in the magnetic moment, constant over a wide 
range of Gd layer thicknesses, was found for the U/Gd system (~ 4/^b compared 
*T , with 7.63/jb for the bulk metal). This could be understood on the basis of a 

, pinning of Gd moments arising from a column-like growth mechanism of the Gd 

layers. A study of the effective anisotropy suggests that perpendicular magnetic 
' anisotropy could occur in multilayers consisting of thick U and thin Gd layers. 

C^") ' A reduction in the Curie temperature was observed as a function of Gd layer 

■ * thickness, consistent with a finite-size scaling behaviour. 

o 

o 

PACS numbers: 61. 12. Ha, 68.35.Ct, 68.65.Ac, 75.30.Cr, 75.30.Gw, 75.70.Cn 
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1. Introduction 

The past two decades have seen a new branch of condensed matter physics develop, 
devoted to the investigation of multilayer thin films, driven principally by the ability 
to manipulate directly the electronic behaviour of materials. The tailoring of the 
multilayer composition and the production of samples consisting predominantly of 
interface regions amplify the interaction of the electronic bands with respect to the 
bulk layers. The study of proximity effects at the interfaces of these multilayer 
systems provides a route into the understanding of fundamental electronic properties. 
Our interest lies in the interaction of the U 5/ electrons with those of the itinerant 
transition metal 3d states of iron and cobalt, and the more tightly bound ^/electrons 
of gadolinium. 

An investigation of the structural and magnetic properties of the U/Fe system 
has been reported previously [TJ [2] , but recent modifications to the growth apparatus 
have allowed the introduction of Nb buffer and capping layers. Paper I in this series 
of articles [3] describes the fabrication and structural characterisation of U/Fe, U/Co 
and U/Gd multilayers grown in this way. 

Earlier work on U/Fe multilayers [2] showed that the magnetism is dominated 
by that of the Fe atoms, and we anticipate a similar situation in the present study. 
Although any polarisation of the U atom is interesting, and is one of the objectives 
of this programme of multilayer growth, we expect any such moment to be at most 
~ 0.2/ib, so that its observation requires the use of element-specific techniques, such as 
X-ray magnetic circular dichroism (XMCD) [4] or X-ray resonant magnetic reflectivity 
(XRMR) [5]. Such techniques are not discussed in the present paper. However, the 
influence of the uranium on the growth of the multilayers, for example its atomic size 
relative to the other atoms in the multilayer structure and the connection between 
growth and anisotropy of the magnetic layers will be reported. As shown in paper I, 
the growth of U/TM (TM = Fe, Co) multilayers results in a poorly defined interface 
which has a large amount of roughness or interdiffusion relative to the thickness of 
the layer. X-ray diffraction spectra indicated that the uranium was present in either 
an amorphous or in a poorly crystalline alpha-U phase. On the other hand, in the 
case of U/Gd both the reflectivity and high-angle X-ray diffraction results show that 
relatively sharp interfaces are produced, and that the uranium adopts an hep form 
with a c-axis dimension not far from that of Gd, which grows in an hep (001) preferred 
orientation. 

This paper will address the bulk magnetic properties and to extend our knowledge 
of the structure within the U/TM and U/Gd multilayers. A combination of 
magnetometry and polarised neutron reflectivity (PNR) techniques have been used to 
investigate the magnetic substructure of the ferromagnetic layers, and for the case of 
the U/Gd system, to study the anisotropy and magnetic ordering (Curie) temperature 
as functions of the bilayer composition. 

2. Bulk Magnetisation 

This section deals with magnetisation measurements carried out on the U/TM 
and U/Gd multilayer systems. SQUID magnetometry measurements probe only 
magnetisation properties from the bulk of the sample, but it is possible to resolve 
effects from each of the individual components by investigating layer thickness, t, 
dependent properties. For the case of the U/Fe and U/Co samples, the dependence of 
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the saturation magnetisation on tp c and tc is presented, while for the U/Gd system, a 
more complete investigation has been carried out. Both ty and t Gd dependent studies 
of the ferromagnetic saturation and of the coercive field are presented. A comparison of 
the field dependent magnetisation measurements taken in the plane and perpendicular 
to the plane of the film gives a measure of the anisotropy present in these multilayers, 
and a zero-field cooled (ZFC) study of the magnetisation as a function of temperature, 
in a constant applied magnetic field, has been used to measure the Curie temperature 
(Tc) as a function of tGd. 

2.1. Experimental method 

The magnetisation measurements were carried out in a Quantum Design MPMS 
(Magnetic Property Measurement System) SQUID magnetometer both at UCL and 
at the Clarendon Laboratory, Oxford. The data were taken using both the DC and 
the reciprocating sample option (RSO) in no-overshoot mode, which enables a slow, 
stabilisation of the field, an important consideration when measuring ferromagnetic 
signals with a large diamagnetic background, as in our case, due to the relatively 
thick sapphire substrates. The measurements of the magnetisation as a function of 
applied field were recorded between -70kOe and 70kOe, at 10K, well below the Curie 
temperatures of iron, cobalt and gadolinium (1043K, 1388K, and 296K respectively). 

2.2. U/TM Field dependence 



Table 1. The saturation magnetisation values as determined by SQUID 
magnetometry for a selection of U/Fe and U/Co multilayer samples. Values are 
included for the absolute magnetisation, scaled by the respective areas and number 
of bilayer repeats and the values of magnetic moment per (Fe,Co) atom. 



Sample Number Composition tF ,Co (±2A) M sat (emu/unit) M sat (/xe/Fe, Co) 



SN71 


[U 9 /Fe 34 ]3o 


34 


1.09E-4 


1.4 ±0.16 


SN72 


[iWFcirJio 


17 


8.22E-6 


0.2 ±0.04 


SN74 


[U 32 /Fe27]30 


27 


7.86E-5 


1.23 ±0.15 


SN75 


[U 35 /Fe27]30 


27 


7.71E-5 


1.2 ±0.15 


SN76 


[U 2 7/Fe57] 2 o 


57 


2.32E-4 


1.72 ±0.19 


SN108 


[U28/CO27]20 


27 


4.38E-5 


0.94 ±0.12 


SN112 


[Ui 9 /Coi 9 ] 2 


19 


9.50E-6 


0.30 ±0.05 


SN114 


[U 9 /Coi 8 ] 2 o 


18 


4.44E-6 


0.19 ±0.05 


SN116 


[Ui 9 /CO42]20 


42 


7.87E-5 


1.10 ±0.23 



Table [T] lists the saturation magnetisation, M sat in units of emu/unit, where the 
values given arc normalised to the respective sample areas and number of bilayer 
repeats. Values of M sat are also given in units of fj,B /Fe, Co. Graphs of these properties 
are plotted as a function of tpe.co hi figure [IJ panel (a) describes the variation of the 
absolute magnetisation, while (b) plots the values of magnetic moment per atom. A 
straight line fit to the data in (a) can be made to describe the linear variation in 
absolute magnetisation, which can then be scaled in the same manner as the raw 
experimental data, to give the trend in saturation moment. At large values of tF e ,Co 
(insert of figure Q] (b)) Ms tends towards the expected bulk moment values for Fe and 
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Figure 1. Variation of the saturation magnetisation, scaled by the area and 
number of bilayer repeats, is presented in panel (a), as a function of the 
ferromagnetic layer thickness. A straight-line fit is shown for both U/Fe and 
U/Co systems. Values of M sa t, given in /ie/Fe, Co, are shown in (b) and a curve 
is presented, calculated from the fit in (a) and extrapolated to large values of 
tFc,Co (insert). Values of M sa t as determined by PNR measurements are also 
presented in (b). 



Co of 2.2/iB and 1.7/l*b respectively. The limited number of data points (figure [T] (a)) 
implies some uncertainty in the extremal values of Ms, although it is clear that values 
close to the bulk would be achieved with very thick layers. Values of the saturation 
magnetisation in /Fe, Co, as determined by polarised neutron reflectivity, are also 
plotted in figure [T] (b); these measurements are described in detail later in the text. 

The linear fit to M sat described in figure Q] (a) shows a direct proportionality 
between the increase in Fe or Co layer thickness and the saturation magnetisation. 
The intercept of this line with the x-axis is an indication of the thickness at which 
there will be no magnetic moment, a magnetic 'dead' layer. However, this distinct 
separation of a non-ferromagnetic component and a ferromagnetic one is an unrealistic 
description of the Fe and Co layers. Mossbauer measurements on the U/Fe system [2] 
have recently been reinterpreted [3] to give an Fe layer comprising of three components, 
Fe amorp hous, which is paramagnetic, Feb cc carrying the bulk moment, and a non- 
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magnetic Fe component likely to be present in a U-Fe alloy, where the spin up and 
spin down 3d bands of the iron can be equally populated; this alloy will be present at 
both U | Fe and Fe | U interfaces, since one of the principal factors for the formation 
of these alloy regions is that of chemical interdiffusion, a process independent of the 
sputtering sequence. 

It is not at first obvious where within the Fe layer these components form, but 
careful consideration of the interfacial structure from X-ray diffraction data [3] and 
results from Mossbauer spectroscopy [2J suggest that the centre of the layers are 
comprised of bulk-like Fe, provided the Fe layers are thicker than the crystalline 
limit, ~ 20 A [3] and then an interdiffuse U-Fe alloy region will be formed at the 
interfaces. These interfaces will then include components of non-magnetic Fe, where 
the concentrations of Fe in U are low, and hence the coordination number of Fe is low. 
Then, across the interface from U to Fe, as the concentration of Fe atoms increases, 
paramagnetic Fe will be found and then ferromagnetic Fe as crystallites of bec Fe form. 
This result is supported by the X-ray diffraction data presented in Paper I [3], which 
shows that for sample SN72, [U23/Fei7]io, there is no visible intensity from any bec 
Fe component, however a small magnetic moment is still present. This result suggests 
that the interface is a complicated mixture of a broad spectrum of Fe components. 
The similar atomic size (10. 3A 3 for Co compared with 11. 5A 3 for Fe) and diffusion 
properties of cobalt are likely to result in a similar interfacial structure and this is 
supported by the results shown in figure Q] 

Although it is not wholly accurate to assign a 'dead' layer value to a multilayer 
system, it is a useful tool to compare the extent of the interface from U/Fe to U/Co 
and to U/Gd. In the case of U/Fe and U/Co the values are very similar, tp° ad and 
tco ad are both ~ 15A. The formation of such interfaces in transition metal multilayers 
is not uncommon and has been reported in Ni/Cu [6], Fe/W [7], Fe/V [8], Fe/Nb [9] 
and Co/Ti [10] systems. 

2.3. U/Gd Field dependence 

A study of the structural properties of U/Gd multilayers [3], compared with those of 
U/TM systems, indicated that the U/Gd interfaces were sharper than those for U/Fe 
or U/Co multilayers and did not include such a significant region of interdiffusion. In 
this case, the magnetic 'dead' layer was expected to be small. 

Values of the saturation magnetisation, M sat , are summarised in table in both 
units of emu/unit and /xe/Gd, where emu/unit refers to a normalisation of M sat per 
unit area and per bilayer. These data are shown in figure [2] (a) for samples with 
tu ~ 25A. A linear relationship is found between the saturation magnetisation and the 
gadolinium layer thickness. The x-axis intercept of the fitted line gives the thickness of 
the magnetic 'dead' layer, which in this case is ~ 4A. Figure[2](b) shows the variation 
of M sat as a function of ted in units of /ib / Gd. The straight-line fit to the data in panel 
(a) has been scaled in the same manner as the experimental data. The inset figure in 
panel (b) shows the extrapolated form of the saturation magnetisation for large values 
of tod, towards a value of 4.3/ib- Figure [2] (b) also includes saturation magnetisation 
data obtained from polarised neutron reflectivity data (PNR), described in Section 3. 

There are some important observations to be made from figure [21 First, the 
magnetic 'dead' layer of about 4A is significantly thinner than that observed in the 
U/TM systems. Second, whereas for the U/TM multilayers M sat tended towards the 
bulk moment value for very thick TM layers (see insert of figure [T] (b)) the saturation 
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Table 2. Summary of the absolute saturation magnetisation (normalised per 
unit area and to the number of bilayer repeats) and M sa t (/ie/Gd), for all U/Gd 
samples measured. 



Sample Number 


t Gd (±2A) 


tu (±2A) 


M sat (emu/unit) 


M sat (/WGd) 


SN63 


33 


26 


3.67E-6 


3 96 ± 48 


SN64 


54 


20 


6.04E-6 


q 07 4- n 44 


SN65 


76 


2G 


8.63E-6 


4.04 ±0.42 


SN66 


20 


39 


1.74E-6 


3.09 ± 0.44 


SN68 


20 


89 


1.46E-6 


2.60 ±0.37 


SN120 


79 


29 


9.32E-6 


4.20 ±0.44 


SN121 


38 


11 


4.23E-6 


3.96 ±0.46 


SN122 


11.4 


11.1 


1.03E-6 


3.21 ±0.73 


SN124(600K) 


24.8 


10.6 


2.33E-6 


3.34 ±0.46 


SN134 


19.8 


10 


2.29E-6 


4.11 ±0.62 


SN135 


18.2 


15.8 


2.02E-6 


3.95 ±0.56 


SN136 


19.4 


19.2 


2.16E-6 


3.96 ±0.56 


SN137 


19.5 


28.2 


1.96E-6 


3.57 ±0.54 


SN138 


20 


4.8 


2.67E-6 


4.75 ±0.67 



magnetisation extrapolated for large values of ted in U/Gd multilayers tends towards 
a value of 4.3^b, significantly reduced from the bulk metal value of 7.63/ib- This 
indicates that the reduction in magnetic moment is not confined simply to the interface 
regions, but is an effect arising from the bulk of the Gd layers. Significantly reduced 
values of the ordered gadolinium moment have been observed in other multilayer 
systems, such as Gd/Mo [11], Gd/V [12] and Gd/Cr [13]. The possible mechanisms 
for this phenomenon were discussed for the Gd/Mo system [IT] in terms of the Gd 
growth and pinning of the Gd moments at the interface. 

In order to probe further the effect of the Gd structure on M sat , a thin sputtered 
film of Gd was grown (sample SN62, ~ 500 A thick) under the same conditions of 
temperature, Ar pressure and sputtering power as the U/Gd multilayers. A value 
of only ~ 5/itB was determined; this indicates that the majority of the magnetic 
moment reduction observed in the multilayers arises from the growth mechanism of the 
gadolinium layers and is not confined to the interface region. Recalling discussions 
presented in paper I of this series [3], for thick Gd layers large roughnesses were 
observed; these were described as step-like, a feature related to a columnar growth 
modulation. It is possible that it is at the boundaries of these column structures 
that the majority of the Gd moments are pinned, although some pinning will also be 
present at the interfaces. 

It is clear, however, that the presence of the uranium layers affects the 
magnetisation of the gadolinium layers. In order to highlight this point figure [3] shows 
the saturation magnetisation values (^e/Gd) for a number of U/Gd samples with a 
constant gadolinium layer thickness of ~ 20A. As the U layers become thicker, so the 
saturation magnetisation decays. 

A likely source of this effect is again a structural one, where variations in tu can 
alter the internal strain within the Gd layers, which can in turn affect the pinning 
of the Gd moments near the U/Gd interface. A similar effect was observed in the 
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Figure 2. The variation of the saturation magnetisation, scaled by the area and 
number of bilayer repeats, is presented in panel (a), as a function of tQ^- A 
straight-line fit is shown. Values of M sa t, given in jig /Gd, are shown in (b) and a 
curve is presented, calculated from the fit in (a) and extrapolated to large values 
of tQd (insert). The starred data points refer to sample SN124, [Uio.6/Gd24.8]20i 
grown at 600K. Values of M aa t as determined by PNR measurements (open 
squares) arc also presented in (b). 



Gd/Mo system [14], but in this case a linear relationship was observed between the 
Mo layer thickness and the reduction in Gd moment. However, the range of Mo layer 
thicknesses discussed was only between 7 and 15A, whereas in our case the U spacer 
layer thickness varies between 5 and 90A. 

2.4- Anisotropy 

The magnetic anisotropy is a measure of the preferred direction of the magnetisation 
[15] . and can be determined directly from bulk magnetisation measurements, where 
the magnetic field is applied both perpendicular to and in the plane of the film. In 
this case a selected number of U/Gd samples were cut so that the measurements 
could be made with the samples in both orientations. The resulting hysteresis loops 
give values for the coercive fields, termed H^, ara and H^ 0lp for the two magnetic field 
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Figure 3. Values of M sat are presented as a function of tu in units of fi^ /Gd 
for samples with t(3<j ~ 20A. A PNR measurement is shown as the open square. 
The solid line is a guide to the eye. 



directions, respectively. A discussion of the anisotropy within these systems is allied 
with a consideration of the values for the coercive field, since the energy involved 
in directing the moments along an easy axis is directly related to the applied field 
required to rotate them. However, it should be noted that the coercive field is also 
dependent on a number of other complicating factors, such as domain size and strains 
within the layers. 

The magnetic anisotropy energy, K e g can be determined from the difference 
in area of the magnetisation loops; this relationship holds for all directions of the 
magnetic easy axis. The sign convention adopted for K e g, denotes that a positive 
K e g indicates a perpendicularly oriented magnetic easy axis. This implies that the 
difference in area of the two loops must be taken as the subtraction of the parallel 
from the perpendicular loop. 

The magnetic anisotropy energy or effective anisotropy is primarily a combination 
of magnetic dipolar (favoring an in-planc alignment) and magnetocrystallinc 
anisotropics, stemming from the long range magnetic dipolar and spin-orbit 
interactions, respectively. These contributions can be resolved into volume, Ky, and 
surface, Kg, contributions to give, 

K eff = K v + 2Ks/t Gd (1) 

Ky is dominated by the magnetic dipolar or shape anisotropy, whereas the lowered 
symmetry at the interface leads to a dominant spin-orbit effect and therefore surface 
anisotropy. A plot of the ferromagnetic layer thickness, t, versus tK c ff generally yields 
a linear function, where the gradient of the line can be used to calculate the volume 
anisotropy term and the y-intercept gives twice the surface term. 

Although this treatment of the anisotropy of multilayer systems has become 
commonplace, we note that there are several important assumptions concerning this 
extraction of the volume and surface contributions from the effective anisotropy. First, 
it is assumed that the anisotropy, localised at the interface region influences the 
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magnetic moments within the bulk of the layer; this is only true if the anisotropy 
is much smaller than the intralayer exchange. Also, the validity of the separation of 
the effective anisotropy into surface and volume terms becomes questionable when the 
layers are very thin and are almost entirely comprised of interface region. In this case, 
Ky is taken as independent of the thickness of the films, but it is possible in multilayer 
systems where the lattice mismatch between the respective species is low, to produce 
strain effects throughout the multilayer that introduce a magnctoclastic contribution 
that changes the magnctocrystallinc anisotropy. 

TableOgives H paia , H pcrp and K e g for selected U/Gd samples. Figure|3](a) shows 
a plot of tcdKeff as a function of the gadolinium layer thickness and as a function of 
the uranium layer thickness. Panel (b) of figure [4] presents the coercive field values as 
a function of ty with the applied magnetic field parallel to the plane of the sample 
and perpendicular. 

Table 3. Summary of the coercive fields, determined from the hysteresis loops, 
taken with the field applied parallel and perpendicular to the plane of the film. 
Values of the magnetic anisotropy energy, K c ff , calculated as the difference in area 
of the hysteresis loops, measured in orthogonal field directions are also listed. 



Sample Composition H para (Oe) H pcrp (Oe) K eff (MJ/m 3 ) 

(±10) (±10) (±0.01) 



SN63 


[U26/Gd33]20 


350 


520 


1.09 


SN64 


[U26/Gd54]20 


430 


510 


1.43 


SN65 


[U2 6 /Gd 76 ] 2 o 


550 


810 


0.98 


SN66 


[U3g/Gd2o]20 


320 


160 


0.68 


SN68 


[Usg/Gd2o]20 


540 


220 


0.63 


SN134 


[Uio/Gdig.8]30 


270 


670 


0.92 


SN135 


[Ui5.8/Gdi8.2]30 


250 


700 


1.06 


SN136 


[Uig.2/Gdig.4]30 


240 


220 


1.12 


SN137 


[U28.5/Gdig.5]30 


260 


730 


0.89 


SN138 


[U4.8/Gd2o]30 


280 


2270 


1.26 



In bulk hep gadolinium a small uniaxial anisotropy acts to align the spins parallel 
to the c-axis between the magnetic ordering (Curie) temperature, ~ 296K, and the 
spin re-orientation transition temperature at ~ 230K. The moments then begin to tilt 
away from the c-axis, reaching 60° at 180K and then move back towards a canting 
angle of 30° at low temperatures [16]. The direction of growth of the gadolinium 
layers in these U/Gd multilayers has been well described [3]; the Gd layers are 
oriented [001], such that the uniaxial anisotropy might act to align the moments 
out of the plane of the film. Figure 0] (a) clearly shows that for all Gd thicknesses 
investigated, a negative effective anisotropy was measured. This indicates an in-plane 
alignment of the magnetic moments. These data points can be fitted to a straight 
line (shown as a dashed black line in figure [4] (a)), which has a negative slope and 
a positive y-intercept. These yield values for Ky and Ks of — 1.2 ± 0.3MJ/m 3 and 
0.28 ± 0.2mJ/m 2 , respectively, which are smaller than those obtained for the U/Fc 
system [17]. The negative slope indicates that the dipolar magnetic anisotropy is 
acting to align the moments within the plane of the film, but the positive surface 
contribution suggests the presence of a magnetocrystallinc anisotropy, which acts to 
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Figure 4. (a) tGdK c ff as a function of tod (solid black squares) and as a 
function of ty (solid red circles). The dashed black line is a linear fit to the 
data, (b) coercive field values as a function of tij with the applied magnetic field 
parallel to the plane of the sample, H^, ara , (left-hand y-axis, solid red circles) and 
perpendicular, iiQ Crp (right-hand y-axis, open blue circles). The solid lines are 
guides to the eye. 



align the spins perpendicularly. The line crosses the x-axis at a thickness of ~ 5A, 
thus a [U25/Gc15]n sample might exhibit perpendicular magnetic anisotropy (PMA). 

The errors on the values of Ky and Ks were relatively large and the data points 
at large Gd thicknesses do not fall on the line. This may be due to roughness effects, 
which are known to dramatically alter the anisotropy in multilayers. These have been 
treated theoretically [TH] and have been observed experimentally in the Co/Pt system 
[19] . A variation from linearity at large values of tod is consistent with the observation 
of increased intcrfacial roughness for thick Gd layers [3]. The effect of varying ty on 
K e g is of a much smaller magnitude. However, for samples with thick U layers K e ff 
tends towards more positive values, favouring an alignment of the magnetic moments 
perpendicular to the plane of the film. 

Figure 2] (b) shows the variation in coercive field for different U layer thicknesses 
for both field directions. A minimum in H^ aia and H^ crp can be observed at tu ~ 20A. 
As tu becomes thinner Hji ara changes very little, whereas H{i olp shows a dramatic 
increase. For thicker U layers the situation is the reverse; H porp changes very little, 
but H paia shows a marked increase. This could be a consequence of the increased 
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strain, a possible explanation for the decrease in M sat seen in figure [3] which would 
lead to a more difficult rotation of domains within the plane of the film. 

2.5. Temperature dependence 

The temperature dependence of the magnetisation has been measured for a number 
of U/Gd multilayer samples with different gadolinium layer thicknesses. Zero-field 
cooled (ZFC) measurements were made between 5K and 375K in an applied field of 
IkOc. In this way it was possible to determine the Curie temperature (Tc) of the 
gadolinium layers as a function of Gd layer thickness. 

It is notoriously difficult to precisely determine the ferromagnetic transition 
(Curie) temperature from magnetisation measurements and most attempts are made 
to describe the paramagnetic phase, using the Curie- Weiss law for localised moments 
above Tq. In our case, with such a small quantity of material and a relatively 
large diamagnetic background from the sapphire substrate, measurements of the 
paramagnetic susceptibility at these temperatures would be very difficult. For the 
purpose of this study Tc was taken as the point at which spontaneous magnetisation 
begins to be observed in the temperature dependence, where the rate of change of the 
susceptibility variation is greatest, i.e. the maximum in the second derivative of the 
susceptibility. 

Figure[5] (a) shows the temperature dependent magnetisation curves for a selection 
of U/Gd multilayers. Panel (b) then plots the Curie temperatures for each of these 
samples as a function of ted ■ The gadolinium layer thickness in this case is presented 
in units of monolayers (ML), determined from X-ray diffraction measurements [3J as 
2.9A. Such a decrease in Tc is well-known in thin film systems and can be described 
by a finite-size scaling relationship. A similar study of the thickness-dependent Curie 
temperature of gadolinium grown on tungsten |20j reported a finite-size scaling effect 
and asserted that an observation of this effect in thin films implies a layer-by-layer 
growth. 

T c {bulk) - T c {t Gd ) x 

Tbp) " 6 ° Gd { ) 

Co is an arbitrary constant, which includes contributions from interlayer coupling 
effects and A = \/v, where v is the three-dimensional Ising critical exponent of the 
correlation length. This treatment of finite-size scaling behaviour well describes 
qualitatively, Ni [2T], Fe [22] . Co [33] and Gd [H] systems studied previously, but 
includes some assumptions about the nature of the system studied. The size of the 
magnetic moment is not taken into account and is assumed to carry an equivalent 
value per atom for different layer thicknesses. This requires a coherent multilayer 
growth with little diffusion at the interfaces. 

A fitted curve, consistent with equation (3) for the finite-size scaling, is shown 
as a solid black line in figurc[5] (b) and yields values for A and Co of -1.56 and 
5.75 respectively. This gives a value for the critical exponent, is, of 0.64, which is 
consistent with a 3D Ising model |25j and is consistent with the relationship observed 
in the Gd/W thin films [20]. These values are quoted within an error of 5%, due to 
uncertainties in Tc and ted- It is also noted that the presence of capping layers on 
thin ferromagnetic films [26] can have an appreciable effect on Tc due to modifications 
in the bandstructure from hybridisation effects at the interface. 

For U/Gd multilayers with the thinnest gadolinium layers, the observed Tc is 
much higher than that expected from the finite-size scaling behaviour. This has 
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Figure 5. The temperature dependent evolution of the magnetisation in an 
applied magnetic field of IkOe is shown in panel (a) for a number of U/Gd 
multilayer samples. The black arrow signifies the direction of increasing Gd 
layer thickness and the starred magenta data points arc for the sample grown 
at an elevated substrate temperature. Panel (b) plots the Curie temperature as 
a function of gadolinium layer thickness. The solid black line is a fitted curve to 
a finite-size scaling relationship. 



also been observed in the Gd/W system [20], [27] and has been related to the cross- 
over from 3D to 2D magnetic behaviour, where theory predicts v = 1.00 in the 2D 
regime. Sample SN124, grown at a temperature of 600K, also exhibits a Tc that 
is higher than that expected. This follows observations of the Gd/W system [20] . 
where the increased Tc values for respective Gd layer thicknesses were attributed to 
the accommodation of misfit dislocations and the presence of large inhomogeneous 
strains, caused by steps and other defects at the interface. Previous results have 
shown that as the substrate temperature is increased the crystallinity is improved, 
but the layer roughness increases [3] . This is consistent with a columnar growth and 
step-like roughness profile, which could produce the effects necessary for an elevated 
Curie temperature. 
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3. Polarised Neutron Reflectivity (PNR) 

Neutrons are an extremely useful tool to simultaneously study the chemical and 
magnetic structure, since they interact with both the nuclei and atomic moments. 
Polarised neutron reflectivity can be used to detect the spin-flip (SF, R H and R 
and non-spin-flip (NSF, R ++ and R~~) scattering channels over a typical Q range of 
0.005 - 0.25A- 1 . In this way, PNR can be used to determine the layer thicknesses 
and roughnesses, the average magnetic moment per atom, the orientation of the 
magnetisation and the distribution of magnetisation in magnetic multilayers. 

3.1. Experimental method 

PNR measurements were carried out on the CRISP rcflectometcr of the ISIS time-of- 
flight neutron source at the Rutherford Appleton Laboratories in Chilton, Oxfordshire. 
The CRISP instrument is situated after the liquid hydrogen moderator at 20K and 
analyses a wavelength band of 0.5-6.5 A. The neutron beam arrives at the experimental 
hutch at an inclination angle of 1.5° and the entrance slit provides a cross-section that 
is 40mm wide and between 0.5 and 6mm high. The wavelength band is defined by a 
disc chopper with an aperture and a nimonic chopper is used for pulse suppression. A 
non-adiabatic spin flipper is used to select the neutron spin direction and polarisation 
relaxation is suppressed by a guide field to the sample position. The sample sits 
approximately 10m from the moderator and reflects neutrons through two further 
slits to the detector, 1.75m away. 

U/Fe and U/Co samples were measured at room temperature and PNR data were 
collected in the on-specular geometry in an applied magnetic field of 4.4kOe; large 
enough to saturate the magnetisation within the plane of the film. U/Gd samples 
were cooled (ZFC) to 10K, well below the their respective Curie temperatures. Since 
the measurements were made at magnetic saturation, only the NSF channels were 
collected. The xPOLLY programme [5H] was used to calculate polarised neutron 
reflectivities for given samples and these were fitted to the experimental data. 

3.2. Results 

The following section presents the experimental data and calculated reflectivities for 
a selection of U/Fe, U/Co and U/Gd samples. The layer thickness and roughness 
parameters were fixed at values obtained from X-ray reflectometry. For the case of the 
U/TM multilayers the large roughnesses accounted for the diffusion at the interfaces, 
responsible for the magnetic dead' layer observed in the bulk magnetisation data. In 
these systems the ferromagnetic layers were separated into regions, which exhibit a 
magnetic moment per atom equivalent to the bulk metal, placed at the centre of the 
layers, and regions with no magnetic moment, at the interfaces. Although it is clear 
that the layers will not have such discrete sections, this technique applied to these 
samples was not sensitive to more complex distributions of the magnetisation. 

It is possible to accurately determine the average magnetisation per atom, by 
closely monitoring the splitting of the two spin channels at the critical edge; the point 
at which the neutron transmits through the multilayer and no longer undergoes total 
external reflection. It is in this region of the reflectivity curve that the splitting of the 
two spin channels is sensitive to the magnetisation of the whole sample. The values 
for the average magnetic moment per atom have been displayed in figures [T] (b) , [5] (b) 
and[3j 




Figure 6. Panels (a) and (c) present the polarised neutron reflectivity 
data measured in the specular geometry at 300K and 4.4kOe for samples 
SN71 — [Ug/Fe34]3o and SN118-([Uio/Co34. 5)20), respectively. Experimental 
data and fitted, calculated curves are shown as black points (curve) for the R++ 
channel and red points (curve) for the R . The asymmetry for these two samples 
is shown in panels (b) and (d); red points represent the data and the solid blue 
line, the fitted calculation. 



Figure [6] (a) shows the normalised non-spin-flip scattering from an example U/Fe 
sample, SN71 ([Ug /Fe34]3o). The solid lines represent the fitted calculated reflectivities 
for both spin channels. Panel (b) displays the asymmetry, determined from the 
difference divided by the sum of the reflected intensities. Figure 0(c) and (d) present 
similar results for the U/Co sample, SN118 ([Uio/Co 3 4.5]2o)- 



Table 4. The thicknesses of the three components of the transition metal layers 
for a selection of U/Fe and U/Co samples, determined from the fitted calculations 
of the PNR. 



Sample Composition txMi (A) txM2 (A) txM3 (A) 



SN71 


[U 9 /Fe 3 4]30 


8.0 


22.5 


3.5 


SN74 


[U 32 /Fe 27 ] 3 o 


8.5 


16.2 


2.5 


SN75 


[U 35 /Fe 27 ] 30 


7.5 


15.0 


4.5 


SN116 


[CO42/Uig]20 


11.0 


24.0 


7.0 


SN117 


[U 9 /C0 51 . 3 ]l5 


8.0 


39.0 


4.3 


SN118 


[Ug/Co 3 4. 5 ] 2 


8.0 


22.5 


4.0 
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A model, describing a multicomponcnt ferromagnetic layer in U/TM mul- 
tilayers has been proposed in an earlier investigation of U/Fc multilayers 
[2J, but a new interpretation of these components was given in paper I [3] 
(UTM a noy|TM amor p hous |TM bu ik|TMU al i oy ). Table H gives the thicknesses for a three 
component ferromagnetic transition metal layer. TM1 represents a combination of the 
UTMaiioy and TM am orphous regions and carries no magnetic moment. It has a den- 
sity that is 5% reduced from that of the bulk metal. The central region of the layer, 
TM2, has equivalent magnetisation and density properties to the bulk. TM3 includes 
the TMUaiioy. The calculated reflectivities and asymmetries only well reproduce the 
experimental data when TM1 is greater than TM3; a trend observed for all U/TM 
samples investigated. This is consistent with the current description of the growth of 
U/TM multilayers [3]. The sum of the thickness of layers TM1 and TM3 give the size 
of the magnetic 'dead' layer, ~ 12A, which is close to the value determined by SQUID 
magnctometry (~ 15A). 




Figure 7. The polarised neutron reflectivity data (panel (a)) and the asymmetry 
(panel (b)), measured in the specular geometry at 10K and in an applied magnetic 
field of 4.4kOe are presented for sample SN63 — [U2s/Gd33]20- Experimental data 
and fitted, calculated curves arc shown as black points (curve) for the R + + channel 
and red points (curve) for the R . The asymmetry data is represented by red 
points and the fitted calculation, by the solid blue line. 

Figure [7| and figure [8] show the experimental data and fitted calculations for the 
NSF scattering of samples SN63 ([U 2 5/Gd3 3 ] 20 ) and SN124 ([Ui . 5 /Gd 2 4.8]2o)- The 
latter was grown at an elevated substrate temperature of 600K. 

The U/Gd samples were modelled by a simple bilayer structure. Figure 2 (b) 
showed that the saturation magnetisation did not vary significantly for a wide range 
of Gd layer thicknesses, indicating a constant distribution of the magnetisation within 
the Gd layers. The reduced values of the magnetic saturation calculated from SQUID 
magnetometry measurements were supported by those obtained from the calculated 
reflectivities, fitted to the PNR data. These values arc listed in table [5] for a selection 
of U/Gd samples. 

The input parameters consist of the real and imaginary parts of the neutron 
scattering length, b (xl(U 5 A), the atomic number density, N (xl0 28 atoms/m 3 ), the 
magnetic moment per atom (/Lis/atom), the angle between the moments and the 
applied field, theta (taken to be 0° in our case, since a magnetic field large enough to 
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Figure 8. The polarised neutron reflectivity data (panel (a)) and the asymmetry 
(panel (b)), measured in the specular geometry at 10K and in an applied magnetic 
field of 4.4kOe are presented for sample SN124 — [Uio.s/Gd24.g]20- Experimental 
data and fitted, calculated curves are shown as black points (curve) for the R++ 
channel and red points (curve) for the R . The asymmetry data is represented 
by red points and the fitted calculation, by the solid blue line. 

Table 5. Summary of the input parameters used to calculate the polarised 
neutron reflectivity data, fitted using the xPOLLY 1281 programme. Values for 
the respective layer thicknesses, densities and magnetic moments are included. 



Sample Layer b (xl(T 5 A) N (xl0 28 atoms/m 3 ) t (A) ^ B /Gd 



SN63 


Gd 


6.500 


2.6 


31.1 


3.9 




U 


8.417 


4.8 


25.0 


0.0 


SN65 


Gd 


6.500 


3.0 


79.0 


3.8 




U 


8.417 


4.6 


26.0 


0.0 


SN67 


Gd 


6.500 


3.0 


63.4 


3.1 




U 


8.417 


4.8 


20.0 


0.0 


SN124 


Gd 


6.500 


2.7 


24.8 


3.9 




U 


8.417 


4.8 


10.0 


0.0 



saturate the gadolinium layers was applied), and the layer thickness, t (A). 

It is clear from figures [7] (a) and (b) that the calculated reflectivities for the NSF 
channels do not well reproduce the data far away from the Bragg peak positions. In 
particular, the R intensity is significantly larger in this region than that calculated. 
This results in a large asymmetry in the calculated case, but very little splitting of the 
NSF channels in the experimental data. The saturation magnetisations, determined 
from the splitting of the channels in the vicinity of the critical edge are consistent 
with SQUID magnetometry measurements. These measurements indicated an even 
distribution of the magnetic moment through the Gd layers. It is therefore likely 
that the discrepancy between calculated and experimental PNR data is a product of 
the mechanism for the reduced magnetic moment. Attempts to introduce a complex 
interfacial structure, as was successful for the U/TM samples, did not result in any 
improvement in the fit for the U/Gd multilayers. 
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4. Conclusions 

The magnetic properties of U/TM (Fe and Co) and U/Gd multilayers have been 
determined, by a combination of SQUID magnetometry and polarised neutron 
reflectivity measurements. These results are consistent with ideas developed in paper 
I [3], concerning the structure and interfacial properties of U/TM and U/Gd systems. 

Magnetisation measurements on U/Fe multilayers grown on sapphire substrates 
with Nb buffer and capping layers, exhibited similar magnetic properties to those 
grown on glass [2], indicating that the predominant effects in this system arise from 
the interfaces of these two elements. Qualitatively similar effects were reported on 
the U/Co system. For large tp c and tc OJ values for the saturation magnetisation close 
to that of the bulk metals were predicted. The presence of a magnetic 'dead' layer 
(~ 15 A) for U/TM multilayers was determined by both SQUID magnetometry and 
PNR. However, it was possible to determine a further degree of complexity concerning 
the TM layers, from the PNR data. The experimental data were well modelled by 
slicing the TM layers into regions with no magnetic moment either side of a central, 
bulk-like region. The thicknesses of these slices were consistent with the structural 
description of the U/TM interfaces developed in paper I [3j- 

A 'missing' moment was also discovered in the U/Gd system, determined by 
SQUID magnetometry measurements and PNR. However, the mechanism for this 
reduction is clearly different from that observed in U/TM multilayers. The 'dead' 
layer reported in this case is < 4A; too small to account for such a large loss in 
magnetisation. In fact, this reduction was constant across a wide range of Gd layer 
thicknesses and was even observed for sputtered thin films of gadolinium. Hence, this 
is not an effect confined to the U/Gd interfaces, but one that is present throughout 
the whole of the Gd layers. A column-like growth, capable of the necessary strains 
and resultant defects in order to pin large numbers of Gd moments is the most likely 
explanation for this reduction; a picture, which is consistent with observations of the 
roughness and crystalline properties described in paper I [3] . 

A study of the anisotropy in U/Gd multilayers was carried out as a function 
of both tu and ted- A precise description of these trends is complex, since a large 
number of effects can alter the magnetic anisotropy; these effects, which commonly 
relate to the thickness of the respective layers, are often correlated to one another. 
However, the dependence of K ff upon the Gd layer thickness, clearly results in a 
negative volume contribution and a small positive surface contribution (although the 
error bars are large), favouring in-plane and out-of-plane alignments of the magnetic 
moments, respectively. The general trend as a function of tu is for an out-of plane 
alignment for thick U layers. Samples with thin Gd layers and thick U present a 
possible route towards PMA. 

A close inspection of the PNR data and calculated reflectivities for the U/Gd 
system revealed a discrepancy between them at the half-Bragg peak positions. This 
feature was represented by an unusually large intensity of the R channel, compared 
with that calculated. Explanations at this point stem from a discussion of the 
mechanisms responsible for the large reduction in saturation moment; that the pinned 
moments at the boundaries of column-like gadolinium grains provide a certain amount 
of diffuse scattering, which could contribute to the intensity observed in the specular. 
Further PNR measurements are planned on thickness matched U/Gd samples, in 
order to suppress the intensity from the even order Bragg peaks and improve the 
sensitivity to the magnetic scattering at these positions. Off-specular measurements 
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are planned to investigate the diffuse scattering so that a model can be developed, 
which can explain such a large reduction in M sat . A programme of measurements is 
also underway to investigate the magnetism of uranium in these multilayer systems, 
utilising the element-selectivity of synchrotron radiation. 
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